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Abstract. It is well known that annihilations in the homogeneous fluid of dark matter (DM) 
can leave imprints in the cosmic microwave background (CMB) anisotropy power spectrum. 
However, the relevance of DM annihilations in halos for cosmological observables is still sub¬ 
ject to debate, with previous works reaching different conclusions on this point. Also, all 
previous studies used a single type of parameterization for the astrophysical reionization, 
and included no astrophysical source for the heating of the intergalactic medium. In this 
work, we revisit these problems. When standard approaches are adopted, we find that the 
ionization fraction does exhibit a very particular (and potentially constraining) pattern, but 
the currently measurable r re i 0 is left almost unchanged: In agreement with the most of the 
previous literature, for plausible halo models we find that the modification of the signal with 
respect to the one coming from annihilations in the smooth background is tiny, below cosmic 
variance within currently allowed parameter space. However, if different and probably more 
realistic treatments of the astrophysical sources of reionization and heating are adopted, a 
more pronounced effect of the DM annihilation in halos is possible. We thus conclude that 
within currently adopted baseline models the impact of the virialised DM structures cannot 
be uncovered by CMB power spectra measurements, but a larger impact is possible if pecu¬ 
liar models are invoked for the redshift evolution of the DM annihilation signal or different 
assumptions are made for the astrophysical contributions. A better understanding (both 
theoretical and observational) of the reionization and temperature history of the universe, 
notably via the 21 cm signal, seems the most promising way for using halo formation as a 
tool in DM searches, improving over the sensitivity of current cosmological probes. 
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1 Introduction 

The existence of a dark component of matter in the universe, i.e. a non electromagnetically 
interacting form of matter, is by now well established thanks to a variety of observation in 
both astrophysics and cosmology: This dark matter (DM) is necessary for instance to explain 
the formation of structures in the universe as we see them, and its relic density Ddm^ 2 can 
be very precisely measured thanks to the cosmic microwave background (CMB) anisotropy 
power spectra. The quest for pinning down DM nature is however still underway, with a wide 
variety of techniques. 

Many extensions of the standard models of particle physics (including electroweak scale 
supersymmetry) naturally accommodate a (quasi-)stable weakly interacting massive particle, 
or WIMP, that can act as excellent dark matter candidate. Additionally, WIMP residual 
annihilations (or, in some models, decays) can inject sufficient “visible” energy that can be 
searched for. Notably, high energy cosmic ray and gamma rays fluxes are routinely analysed to 
reveal excesses that could be attributed to a DM origin, or conservatively to constrain particle 
physics parameters such as the annihilation cross-section times velocity (av) (averaged over 
the velocity distribution, hereafter called “the cross-section”). 

Interestingly, it has been realized since the eighties that CMB observations can in fact 
tell us more about the nature of DM. First of all, purely gravitational arguments can lead 
to robust constraints on its lifetime, independently of the particle physics models (see e.g. 
[1-4]). But for more specific classes of DM candidates, such as WIMPs, the CMB diagnostic 
power is much stronger. Annihilations (or decays, neglected from now on) inject non-thermal 
photons and electrons in the intergalactic medium (IGM) that can delay the recombination 
and change the relic abundance of free electrons after decoupling. Hence, WIMP annihilations 
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can jeopardize the observed CMB temperature and polarization anisotropy angular power 
spectra and therefore can be constrained by an experiment like Planck. 

DM annihilations in the homogenous smooth background have been well studied and 
documented in the last decade [5-23]. The most realistic calculations for WIMPs have been 
done by Refs. [13, 15] and updated recently in Ref. [24], where authors carefully computed 
how much of the initial DM particle energy is deposited into the medium, as well as how 
this energy is separated between ionization of hydrogen atoms, excitation of these atoms and 
heating of the plasma. They also found that the impact of DM annihilations depend sizably 
on the mass and the produced particles (electrons, quarks, etc.). The Planck collaboration 
in a very recent paper [25] has reported very strong bounds on the cross-section, excluding 
thermal WIMPs for any standard model annihilation channel for masses up to 10 GeV, also 
ruling out WIMP explanations of cosmic ray lepton spectral features discussed in recent years. 

Previous CMB studies mostly focused on the impact of annihilations in the averaged 
cosmological density held of DM. However, at relatively low redshift, the DM fluid clusters 
under the action of gravity into virialised structures, so-called “DM halos”. This process 
increases the averaged density square ( p 2 ) with respect to the square of the smooth background 
density, ( p ) 2 , while the two are nearly equal at high redshift. One could naively expect that 
this results in a large enhancement of the annihilation rate and therefore in a significantly 
bigger impact of DM annihilation on the CMB power spectra. But the effects of halos are more 
subtle since, as we will see, the way in which energy is deposited into the medium changes 
as well. Thus, the modification of the bounds on DM annihilation cross-section cannot be 
trivially obtained. 

Another interesting feature of DM halos is their possible impact on the ionization his¬ 
tory. In the standard picture, it is assumed that stars are the only reionization sources. 
Unfortunately, our knowledge of first stars formation in the universe is very rudimentary, and 
hence also is our knowledge of the ionization history. The formation of halos, if it increases 
significantly the DM annihilation rate, could introduce a new source of reionization in the 
universe and leave a very peculiar imprint on the history of the ionization fraction x e (z) and 
temperature of the IGM Tm, also referred to as the matter temperature. In the past, this has 
been invoked as a way to solve a tension on the measurement of the reionization optical depth 
T r eio coming from WMAP data [26] (preferring a high value of r re io and therefore a relatively 
high ionization fraction at redshift z > 10) and the so-called Gunn-Peterson effect as it is 
measured in astrophysics [27-29] (requiring a relatively high neutral hydrogen fraction above 
z — 6.5 and hence pointing towards smaller values of r re i Q ) [10, 11]. The new Planck data in 
the most conservative case yield r re i 0 = 0.079 ±0.017 (Planck TT,TE,EE±lowP, 95%CL [25]) 
and therefore have reduced this tension to a ~2 a level, since now a single-step reionization 
ending at z = 6.5 is marginally compatible with Gunn-Peterson bounds [25]. Yet, it remains 
interesting to quantify the potential contribution of DM halos to this observable. 

Only a handful of articles have investigated the impact of DM annihilation in halos, 
notably [10-12, 14, 20, 22, 30, 31]. Unfortunately, previous authors follow different formalisms 
and are difficult to compare with each other. More importantly, they arrived at different 
conclusions. For instance, the impact of annihilation on reionization is substantial even 
for baseline parameters according to Ref. [20], relevant for light particles models according 
to [14, 30]; while Refs. [10],[11], [12], and [22] find it to be negligible. One of the few points on 
which all agree is that DM annihilation in halos cannot be the only source of reionization in the 
universe: even in ref. [20], an astrophysical contribution is needed at least to account for the 
Gunn-Peterson observations. What is less well understood is the role that DM annihilation 
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in halos can play in mixed reionization scenarios, perhaps easing tensions between Gunn- 
Peterson observations and CMB ones. In order to clarify this situation, we first evaluate 
the impact of DM annihilation in halos onto the ionization history and compute the optical 
depth to reionization T re i 0 in models with conventional reionization from stars: We find that 
halos can play only a minor role, confirming the conclusions found in most previous literature. 
We check that Planck data do not alter the sensitivity to the halo contribution. To clearly 
explain why it is so, we also characterize the impact of the annihilation in halos onto the 
CMB temperature and polarization power spectra, thus updating the study [14]. In doing 
that, we identify and correct a few mistakes and oversimplifying assumptions used in the 
previous literature [14, 20]. We also provide two major improvements over previous works: 
first, by adopting a phenomenological model for the star formation rate and corresponding 
injection of high-energy photons, we study the dependence of the signals of interest from the 
astrophysical model adopted. In fact, till now all works have studied the problem within a 
single type of parameterization for the astrophysical reionization history. Second, we amend 
an unjustified simplification in the treatment of the IGM temperature evolution, by adding an 
astrophysical source term reflecting the corresponding one in the reionization history. Finally, 
we discuss how and why both the ionization history and IGM temperature evolution provide 
more promising perpectives as DM probes via the halo term. 

Our paper is structured as follow. In section 2, we present our formalism. We review the 
standard recombination and reionization equations, as well as the computation of the standard 
energy deposition functions f(z) necessary to quantify how much energy is deposited into the 
medium, and describe how to compute them in presence of halo formation. Subsection 2.1 
is particularly original, and indulges way more than previous literature on several important 
aspects of the astrophysical source terms for the equations of x e (z) and Tm(z). In section 3, 
we compute the impact of annihilations in halos on the reionization history and revisit the 
question about the possible contribution of halos to solve the slight tension between CMB 
and Gunn-Peterson concerning the reionization optical depth r re i 0 . In section 4, we present 
our results concerning the impact of DM annihilation in halos on the CMB power spectra. 
Section 5 contains a summary of our results, as well as a discussion of possible observables 
where DM halos can play a non-negligible role, that would be worth studying in the future. 
For the sake of a self-contained treatment, more details on the standard Peebles equations are 
available in appendix A. Some remarks on the energy deposition functions and a complete 
comparison between the formalisms followed by different authors is developped in appendix 
B. Appendix C summarizes our treatment of halo formation. Finally, appendix D contains 
a discussion of the reionization optical depth r re i 0 measured by Planck, compared to the real 
optical depth to reionization, which aims at justifying (to the best of our knowledge, for the 
first time) within which errors one can assume them to coincide. 

2 Ionization and thermal evolution equations 

Throughout this study we work with the Boltzmann code CLASS 1 [32], in which the two 
recombination codes RECFAST [33] and HYREC [34] are implemented. RECFAST uses fitting 
functions to take into account the effect of highy-excited states. To deal with DM halos, 
RECFAST cannot be used without extensive modifications, for its fitting functions would 
need to be extrapolated beyond their range of validity. We therefore choose to use only 
HYREC, in which those effects are computed at a more fundamental level without requiring 
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interpolation. While accounting for highy-excited state corrections, HYREC (like RECFAST) 
still contains a system of coupled differential equations with the same basic structure as in 
Peebles case B recombination model. For the evolution of the hydrogen ionization fraction 2 
x e , the code integrates 


dx e (z) 

dz 


1 

(1 + z)H(z) 


( R(z)-I(z )) , 


( 2 . 1 ) 


where R and I are the recombination and ionization rates, and for the evolution of the matter 
temperature, 


dTyi 

dz 


1 

1 + 2 


2Tm + 7(7m - Tcmb) 


( 2 . 2 ) 


Details on these terms and coefficients are given in Appendix A. 


Energetic particles injected by any type of sources will have three effects on the cosmic 
gas: direct ionization, collisional excitation (followed by photoionization by CMB photons), 
and heating. These effects are taken into account by adding two terms to equation (2.1) 
and one term to (2.2). The rate of the first two effects, namely direct ionization Ixi and 
excitation+ionization Ixa- are given by: 


Xi(z) dE 
nu{z)Ei dVdt 


(1 -C)x a (z) dE 
nu{z)E a dV dt 


dep 


(2.3) 


where and E a are respectively the average ionization energy per baryon, and the Lyman-a 
energy (see Appendix A for the definition of the quantity C). In practice, we simply define 
an effective ionization rate Ix(z ) = Ixi(z) + Ixa{z ) and add it to I(z) in equations (2.1). 
The heating rate term K b . normalized to the Hubble rate, to be added in the square brackets 
at the RHS of Eq. (2.2), can be similarly defined as: 


R = _ 2Xh(z) _ dE_ 

h H(z)3k b n H (z)(l + fm + x e ) dV dt 

In previous equations, g^jldep stands for the energy deposited in the plasma at redshift 
2 . It is splitted according to the energy repartition fractions Xj , with the index j = {i, a, h} 
denoting ionization, excitation (through Lyman-a transition) and heating, respectively. In 
general, the factorization into a “universal” deposited energy factor times the coefficients 
X’s (only functions of redshift) is not exact, depending on the type and energy of particles 
responsible for the heating, excitation and ionization. However, it is a useful approximation 
for a specific source term (weak scale DM, stellar sources, etc.) and we will adopt it in the 
following. 


(2.4) 

dep 


2.1 Astrophysical source terms 

In absence of DM, the high-z evolution of x e and Tm is the result of solving Eqs. (2.1), (2.2), 
without additional source terms. However, the resulting evolution would be clearly unphysical 
in the range 2 < 0(10). At very least, we know that the low -2 universe is ionized and relatively 
hot. In cosmological applications, (the bulk of) this is implicitly attributed to unspecified 

2 The code also follows the Helium ionization fraction, but like in most of the literature we neglect here 
the impact of DM annihilation on Helium recombination: it has been checked explicitely in [20, 21] that this 
effect is negligible. 


4 
















astrophysical sources, either unaccounted for or described by some prescription by hand. 
For instance, in the paper [20] a prescription similar to the default “single-step reionization” 
incorporated in the CLASS and CAMB codes is used: below some arbitrarily chosen value 
z re io• %e( z ) is cut and matched continuously to a half-hyperbolic tangent centered on z re i 0 , 
reaching an asymptotic value of one (if He is neglected) at z = 0, with a narrow width 
parameter (6z = 0.2) to describe a fast reionization. Apart from using such a modified x e (z), 
no modification at all is included in the evolution of Tm, i.e. no astrophysical sources of 
heating are considered. Basically all previous treatments have followed a similar approach. 
Here we attempt for the first time to quantify the effect of these approximations, comparing 
them with a more realistic treatment of the astrophysical source terms, inspired by recent 
literature. 

It is usually accepted that a dominant source of reionisation would be given by Lyman 
continuum photons from UV sources in pristine star-forming galaxies. To account for these 
photons, we add a source term taken from [35, 36] of the form of Eq. (2.3) to the evolution 
equation of x e : 

= A* /esc£ionPSFR,( 2 )(l + z) (2-5) 

dep 

where f esc is the fraction of photons produced by stellar populations that escape to ionize the 
IGM, £i on is the Lyman continuum photon production efficiency of the stellar population and 
PSFR is the comoving star formation rate density. The fiducial values for these parameters are 
taken from Ref. [36], namely f esc = 0.2, £; on = 10 53 ' 14 s _ 1 M)7, 1 yr. We also use the functional 
form for the star formation density rate provided in Ref. [36], 


1 dE 
EdVdt 


PS FR = 


(1 ± z) b P 

% + [(! + z)/c p ] d v 


( 2 . 6 ) 


with fit parameters a p = 0.01376±0.001 M 0 yr _ 1 Mpc~ 3 , b p = 3.26±0.21, c p = 2.59±0.14 and 
d p = 5.68±0.19 [36]. Rather than varying the fitted parameters within their confidence level, 
for simplicity we leave the overall fudge factor A* (expected to be of order one) free to match 
the measured optical depth. A consistent treatment of such a source term in an equation 
analogous to Eq. (2.3) requires to specify the corresponding energy repartition functions, %’s. 
We follow the prescription of Refs. [5, 20] in which Xi = Xa = (1 — x e )/3. To motivate this 
choice, we remind the reader that in the simple estimate provided by Shull and Van Steenberg 
[37] approximatively 1/3 of the energy is effectively available for ionization in a neutral gas 
{x e = 0). The adopted expression also fulfills the obvious physical criterion that no energy 
is available to ionise (or excite) an already fully ionised gas (x e = 1). A linear interpolation 
is used for values in between, which corresponds to the reasonable Ansatz that the rate of 
ionisation is proportional to the abundance of neutral hydrogen. The equality Xi = Xa has 
been checked to be approximately true by Ref. [5]. 

Qualitatively, stellar phenomena should also contribute to the heating of the IGM. In 
principle, one might expect this phenomenon to be capture by a term similar to the one 
introduced above, modulo a different energy repartition fraction function. However, it is 
sometimes argued that heating is most efficiently achieved thanks to harder photons (X-ray 
band), see e.g. Ref. [38]. To account for this, we introduce a source term in the brackets 
at the RHS of Eq. (2.2) of the form of Eq. (2.4), with a different normalization, taken from 
Ref. [38], 

= 3.4 x 10 10 /a-Psfr(z)(1 + z) 3 ergs _ 1 Mg 1 yr. 
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dVdt 


dep 


(2.7) 







The X-ray efficiency fudge factor fx- expected to be of order 0(0.1), is set to the benchmark 
value fx = 0.2. Concerning the heating repartition function, for consistency with the approx¬ 
imation used above for Xi and Xa we adopt Xh = (1 + 2x e )/3 as suggested by [5, 20]. To our 
knowledge, it is the first time that a term like this is accounted for in a study of cosmological 
sensitivity to DM effects. 

Finally, note that similar modifications could also be done to the equations describ¬ 
ing helium reionisation, but we leave that for future works since HYREC cannot be used in 
its current form to model such subleading effects. For the present study, we keep using a 
phenomenological hyperbolic tangent function to describe helium reionisation. 


2.2 Dark Matter annihilation in the smooth background 

Concerning the energy repartition functions from DM, in this work we use the recent calcula¬ 
tions of Xj ( z ) by Galli et al. [21], that improved over former estimates [5, 37]. An alternative 
computation has been done in Evoli et al. [39], but that of Galli et al. [21] more closely 
resembles the formalism adopted here. The energy density injection rate gf^linj can be read¬ 
ily computed as the product among the number density of pairs of DM particles n pa i rs , the 
annihilation probability per time unit P ann , and the released energy per annihiliation E ann : 


cLE 


dVdt 


(0 — ^flpairs — K 2 ^ ' ^ann — (o'ann'c)^'DM^ ' ^F/ ann — 2 ??IdmC ^ • (2.8) 


Taking only into account the smooth cosmological DM distribution, we can write this rate as 


dE 

dVdt 


(0 

inj,smooth 


KPcC 2 Q dm (1 + ~) 6 


(<7 an nf) 

moM 


(2.9) 


In the equations above, (a ann v) is the cross-section, bdm = Pc^dm(1 + £) 3 the number density 
of DM particles, p c = 3 Hq/8itG the critical density of the Universe today, the current 
DM abundance relative to the critical density and ttidm the DM mass. If DM is made of 
self-conjugated particles, such as Majorana fermions, one has k = 1, which is what we shall 
assume in the following; if DM particles and antiparticles differ (as in the case of Dirac 
fermions) and are equally populated, k = 1/2 since only half of the pairs that one can form 
(the ones made by one particle and one antiparticle) are suitable for annihilation. 

The response of the medium to energy injection depends strongly on the cascade of particles 
produced by DM annihilation, and on the epoch at which the DM particles annihilate. This 
response is conveniently parametrized by a dimensionless efficiency function f(z) [13] such 
that: 


dE 

dVdt 


CO = /(0 


dE 

dVdt 


(0 

inj, smooth 


( 2 . 10 ) 


The expression of / can be obtained via appropriate transfer functions T^\z', z, E), giving 
the fraction of the ^—particle’s energy E injected at z' that is absorbed at z, as 


/(0 


f dln(l + z') 


(l+z') 3 

H(z') 


ZjTW(z',z,E)E^\W d E 


(i+U 3 

H(z) 




ln J 


( 2 . 11 ) 


where the sum runs over species (in practice £ denotes either photons or electrons), and |fc 
is the injected spectrum of each of them in a given DM model, and is independent from z. 
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The calculation of these functions is very involved, but it has been carefully done in [13, 19], 
with a study of associated systematics presented in [21]. While we do not indulge here in 
technical details, it is worth stressing a few conceptual issues concerning the meaning of f(z): 

• In the literature, the assumption that the energy released by annihilations at a given 
redshift is absorbed at the same redshift is referred to as the on-the-spot approximation. 
In that case, the meaning of f(z ) is clear: it is the fraction of energy that is absorbed 
by the gas, either via collisional heating or atomic excitations and ionizations. It takes 
into account that a part of the energy may escape for instance in the form of neutrinos, 
to some extent protons [40], or as photons which free-stream to the present day. The 
f(z ) factor, in this approximation, depends on DM particle model but cannot exceed 
1 by definition. This approximation therefore mainly consists in considering that all 
absorption processes are very rapid in comparison to the Hubble time, defined as tn(z) = 
c/H(z). 

• However, the authors of [13] have shown that this is not true for the entire redshift range 
that we are considering (and strictly speaking, not even at z ~ 1000). Some photons that 
are free-streaming at some given redshift z' could be absorbed at z < z' . The beyond- 
on-the-spot treatment consists in computing the full evolution, like in Refs. [13, 19, 21]. 
The result can still be cast in the form of an efficiency function f(z), simply defined 
a posteriori as the ratio of deposited energy to injected energy at the same redshift. 
In Appendix B, we quickly summarise the principles of the calculation performed in 
Refs. [13, 19, 21], and comment on differences with the approach of [10, 11, 14, 30, 31]. 

Several authors have shown that the redshift-dependence of f(z) is of very little relevance for 
CMB constraints [16-18, 20]. This is because the main impact of smooth DM annihilation on 
the CMB is to inhibit recombination, enforcing x e to freeze out near redshift z ~ 600 at larger 
values than in standard ACDM. Thus, the effects of DM annihilation is usually parameterized 
by a single quantity p a nn defined as: 

Pann = , where f eS = f(z = 600). (2.12) 

?bdm 

However, it is important to keep in mind that, to correctly describe how halos of DM influenced 
the ionization and thermal history of the Universe, one cannot captured DM effects with 
only one model-independent parameter. Furthermore, in order to be able to use the same 
parametrization, it is necessary to recompute these f(z) functions in the presence of halo 
formation, and this will be the main focus of the next section. 


2.3 Dark Matter annihilation in halos 

The spatial average of the annihilation rate is proportional to the average square dark matter 
density. The main impact of structure formation is to increase this average with respect to the 
smooth background case, by an amount usually parametrized through a boost factor B(z): 

(p 2 )(z) = (l + B(z)) (p) 2 (z). (2.13) 


One now has: 

dE 


2 2 
= PcC 


^dm(1 + z ) 


6 ( a a 


(1 + B(z)). 


dVdt 


inj, smooth+halos 


^DM 


(2.14) 






Several ways to compute B(z) have been proposed. We summarize our approach in Ap¬ 
pendix C. The two key (unknown) physical quantities are the maximal overall boost factor 
due to halos and the epoch for the onset of formation of virialised objects. The simplest 
choice adopted in our model was to choose as free parameters the characteristic redshift Zh, 
related to the time of halo formation (occurring near 2 = 2 Zh), and a parameter /), related to 
the amplitude of the boost factor today (since B(z = 0) = f^e rfc(l/(l + +), see equation 
(C.5)). The range of values explored relies on results found in the literature, see Appendix C 
for quantitative details. The evolution of ( p 2 )(z ) for different values of these parameters is 
shown in figure 1. 



Figure 1. DM squared density vs. redshift 2 for several models of halo evolution. 


Some treatments of the effect of DM (including halos) has been presented in the past: for 
instance, Giesen et al. [20] performed this calculation on the basis of a simplified formalism 
developped by Natarayan [14, 30, 31], accounting only for energy deposition through the 
Inverse Compton Scattering (ICS) effect, and taking as a source for the ICS the energy 
injection function of the smooth case. Here however we adopt a treatment based on a more 
straightforward generalization of the equation (2.10) which is equivalent to the one reported 
in [22], We define f(z ) as 


dE 

dVdt 


dep, smooth+halos 


m 


dE 

dVdt 


inj, smooth+halos 


where now equation (2.16) generalizes as 


/(*) 


fdln(l + z') 


(l+z'f 

H(z') 


(1 + B(z / ))ZJ T W (z',z,E)E^dE 


(1+A 3 

H(z) 




(2.15) 


(2.16) 


It is clear that when setting B = 0 one recovers the standard expressions for the smooth 
contribution. 

In order to assess the impact of the improved calculation of the boost factor and of 
the new Planck data, we computed the energy deposition function for two baseline models 
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with annihilation channel XX e + e - and XX ~ as we H as f° r the two masses 
ttidm = 1,1000 GeV, following the generalization of the method of [19] described above and 
making use of the numerical tools provided by the authors 3 . 

The generalization of the calculation for more DM models would only require to adapt 
the injected spectra. Since we aim at being model-independent, we selected (and limit our¬ 
selves to) these two final state examples for two reasons: first, CMB bounds are particularly 
interesting for them, since light leptonic final states are the most difficult models to constrain 
through other methods; second, they represent two extreme cases for the corresponding values 
of f(z) (high for e, low for /r), and hence they are sufficient to bracket typical constraints, if 
re-expressed in terms of (cx a nn)- 

Figure 2 shows our result for the total f(z ) for each baseline models, with ?77 dm = 1 GeV 
or 1 TeV, and with halo parameters [z^ = 30, B = 10 6 ] or [z^ = 20, £> = 10 12 ], compared to 
the functions computed from annihilation in the smooth background only. Note that at high 
redshift, when B <C 1, our result is asymptotically equal to that obtained in [13, 19] in the 
absence of DM halos. Note also that we only performed the calculation down to redshift 
z = 10, since Ref. [21] does not provide transfer functions below this redshift. Physically, 
with our definition, we do not expect big changes of the deposition function at low redshift. 
Hence we assume that f(z) remains constant below 2 = 10, as shown in Figure 2. If this 
assumption turned out to be inaccurate, our final results would not be much affected, because 
observable effects at low redshift are given by the product of f(z) by the factor (1 + B(z)) on 
which there is a huge uncertainty, and that we treat as a free parameter (see eqs . 2.14 and 
2.15). 




Figure 2. Energy deposition functions in the two baseline models, for two values of the dark matter 
mass and several halo parameters, compared to the case with annihilation in the smooth background 
only. Below z = 10, we assume that these functions remain constant. 


3 Impact of high redshift sources on the reionization history 

There is no compelling reason to invoke an extra ingredient such as DM annihilation to explain 
the reionization history of the universe. Actually, our knowledge of pristine star formation 

, 3 http: //nebel. rc. f as .harvard.edu/epsilon/ 
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and its impact on the ionization history x e (z) is so rudimentary that currently, we can treat 
the ionization history x e {z) caused by star formation almost as a free function, and some room 
for an exotic source of reionization is definitely possible. To illustrate this point, in the left 
panel of Fig. 3 we show two possible reionization histories of astrophysical origin: the green 
curve represents the standard step-like model “put by hand”, while the red curve represents 
a model inspired by actual astrophysical data, as described in Sec. 2.1, and normalized (via 
the parameter A* ~ 3) so that the optical depths for the two models are the same. As far 
as cosmological observations are concerned, they are essentially indistinguishable, as we will 
stress again in the following. The points report constraints from [27-29]. In the right panel 
of Fig. 3 we report the corresponding gas temperature evolution, compared with the CMB 
temperature evolution (purple curve): the blue curve represents the typical approximation 
in which this quantity has been evolved in past literature, with only the feedback for the 
x e evolution accounted for (no heating source term). The green and red curves represent 
the evolution of the temperature if a source term similar to the corresponding one adopted 
for x e is included (green: “sudden” heating, put by hand; red: redshift evolution inspired 
by an actual astrophysical model, see Sec. 2.1). The yellow band represents some indicative 
constraints from ref. [41]. Our aim here is not to determine a viable heating history, rather to 
show the rudimentary status of these treatments (with large uncertainties in the astrophysical 
term) and the large room for exotic sources of heating. 




Figure 3. Evolution of x e (z ) (left panel) and Tm(z ) (right panel) in the different approximations 
described in the text, for two prescriptions for describing the effect of astrophysical sources. 


Despite the somewhat unsatisfactory situation, some consensus has been reached on 
important points concerning the reionization history. For instance, in the past the question 
has been raised if the totality of the reionization related phenomenology could be accounted 
for by DM only, but it is now acknowledged to have a negative answer. Even in Ref. [20], which 
finds potentially large effects at high redshift due to DM in halos, an astrophysical contribution 
is needed to account for the Gunn-Peterson effect, requiring the presence of a non-negligible 
neutral hydrogen fraction at redshift z ~ 6.5. On the other hand, CMB observations need 
the Universe to be significantly ionised at higher redshift, in order to get a correct integrated 
optical depth to reionization r re i 0 , compatible with measurements of the temperature and 
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polarization spectra 4 . Although early measurements (notably by WMAP 1 [42]) hinted to 
the necessity for non-trivial reionisation history (e.g. multiple stellar populations or exotic 
DM contribution) because of a tension between these two observables, this is by now mostly 
gone: latest Planck data [25] prefer a lower value of r re i 0 and even a single-step reionization 
ending at z ~ 6.5 is marginally compatible with Gunn-Peterson bounds, reducing the tension 
to a 2 -ct effect. In a single step model, the measurement r re i 0 = 0.079 ± 0.017 (Planck 
TT,TE,EE+lowP, 95%CL [25]) translates for instance into z re i 0 = lO.O^J'g. 

However, mixed reionisation scenarios involving a relevant DM role are still of interest: 
for instance, one may wonder to what extent CMB upper bounds on r re i 0 may lead to stronger 
constraints on DM annihilation than those coming from the smooth background. This was 
for instance the conclusion found in Ref. [20]; even in Ref. [14], it was argued that at least 
light (few GeV’s) thermal relics may have measurable effects. Articles such as [10], [11], 
[12], [22] find it to be instead negligible, for comparable choices of parameters. We want to 
reconsider this with a state-of-the-art approach, correcting some errors and going beyond the 
approximations that we have identified in Refs. [10, 11, 14, 20] as mentioned in Sec. 2.2 and 
developped at the end of Appendix B. We also want to compute the full CMB power spectrum, 
while Refs. [10, 11] only estimated r re i 0 . This extra step may be instructive in establishing to 
what extent future data may improve over current constraints, a possibility raised for instance 
in Ref. [14]. In the following, we fix the key parameter p arm close to the current 95%CL upper 
bound inferred from Planck TT,TE,EE+lowP data 5 [25] with annihilation in the smooth 
background only, namely, p ann = 2.3 x lCU'm 3 * /s/kg= 4.1 x lCD 28 cni 3 /s/[GeV/c 2 ]. If for this 
maximal value, we hud the role of DM annihilation in halos to be negligible, then it will be 
a fortiori true for any viable p ann value. 

As a first sanity check that DM annihilation cannot have a dominant role in a consistent 
reionization history, having fixed p ann to its maximum value, we vary the halo formation 
redshift zp and find the value of the boost factor B{z = 0) giving a reionization optical depth 
of r re io £ [0.045,0.0113] (the 95%CL interval inferred from Planck TT,TE,EE+lowP data). 
Even if the values of the parameters depends on the DM mass and annihilation channels, we 
find that reaching the minimal allowed r re i 0 — 0.079 while assuming the maximal p a nn requires 
halos to form very early and to be very concentrated, e.g., [zp = 50, B[z = 0)) = 10 11 ] for 
a 1000 GeV DM annihilating into muons or [zp = 40, B(z = 0)) = 10 10 ] for a 1 GeV DM 
annihilating into electrons. Cosmic ray data [46] and N-body simulations [47, 48] are hardly 
compatible with B(z = 0) > 10 8 . Hence, reionization from DM annihilation would require 
even greater halo formation redshifts than the maximal value we consider: zp > 50. Even if 
very little is known about the first halos in the universe from the observational point of view, 
such early halo formation times do not appear realistic and are in general not considered 
in the literature. We can conclude in agreement with previous studies that for conventional 
assumptions on annihilating DM models and on halo formation, DM annihilation cannot play 
a dominant role in reionizing the Universe and can at most coexist with stellar reionization. 

A posteriori, this justifies our choice to fix the value of p ann : If plausible reionization 
models involving DM annihilation could produce an exceedingly high r re i 0 , we could use the 
measurement of the optical depth by the Planck satellite to derive new upper bounds on 

4 Strictly speaking, the parameter called reionization optical depth by the Planck collaboration is the 
physical optical depth to reionization only within some assumptions (see appendix D for details). 

5 We checked that the result is exactly the same with the analysis pipeline used in the Planck paper [25], 

based on camb [43] and cosmoMC [44], or using the class version used as a baseline in this work (version 2.4.3) 

[32] and MontePython [45]. 
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Figure 4. Ionization fraction x e (z) as a function of redshift for several mixed reionization models. 
Green lines are benchmark, purely astrophysical reionization scenarios (solid: single step; dot-dashed, 
phenomenological). The blue and red versions of the corresponding lines show the case where both 
smooth DM injection and halo one have been added as well, with growing role of halos, respectively. 
We assume p ann fixed to its most optimistic value. 


this parameter. Since this is not the case, and since we have no precise information on star 
formation, we can always obtain the correct r re i 0 by assuming the maximal realistic effect 
from DM annihilation and a complementary effect from stars. The only hope to obtain new 
bounds on DM annihilation from CMB observations is to analyse the full shape of the CMB 
temperature and polarization spectra: this will be the topic of section 4. 

It is interesting to explore a bit further mixed reionization scenarios. In Fig. 4, we 
show in green the benchmark, purely astrophysical reionization scenarios (solid: single step 
with z T eio = 6.5, as suggested by the Gunn-Peterson bound; dot-dashed, phenomenological). 
The blue and red versions of the corresponding lines show the case where both smooth DM 
injection and halo one have been added as well, with growing role of halos, respectively. We 
have fixed there p ann to the maximum value allowed by Planck [25]. The free parameters are 
the halo ones, z h and 13(0), besides the DM model (annihilation channel and mass), here fixed 
to the muon final state and 1 GeV mass. 

In Figure 5, we vary these two categories of parameters, and find r re i 0 as a function of 
them (in the step-like reionization scenario), together with 68% and 95% confidence limits 
from Planck TT,TE,EE+lowP data, and bounds on B(z = 0) inferred from N-body simu¬ 
lations by Ref. [47]. Figure 5 shows that in absence of DM annihilation in halos, there is 
a marginal (~ 2cr) tension between the model without DM annihilation in halos and the 
Planck bounds on r re i 0 . At the same time, we can also see that DM annihilations in e + e - 
and for realistic values of the boost factor (blue band) do not significantly enhance 
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r re i 0 . The conclusions would be similar for other annihilation channels (all being bracketed 
by these two) or masses. Note that in all these models, we decided to saturate the CMB 
bound coming from annihilation in the smooth background. This means that we fixed the 
observable quantity p ann , but not the fundamental parameter {crv}. This explains why in 
Figure 5, the effect seems to be stronger in the muon case than in the electron case; if the 
cross-section were fixed, the conclusion would be opposite. In this observable, the result is 
also independent of the reionization model adopted, provided that they are responsible for 
the same optical depth. In all cases, our main conclusion is that one needs to push the halo 
contribution to the same unrealistic values as before to remove the marginal tension between 
CMB and the Gunn-Peterson bound. 




Optical depth for several DM models, <^=30, z reio =6.5 
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Figure 5. Reionization optical depth T re i 0 in mixed reionization models, for different DM masses, 
annihilation channels and halo parameters. p ann has been fixed to its most optimistic value and the 
redshift od reionization from stars to = 6.5. The red stripe shows the most conservative bound 
from Planck on r re i 0 [25] and the blue one the most conservative interval for the values of the halo 
amplitudes at z = 0 according to Refs. [47, 48]. 


In summary, our main conclusion is that considering DM annihilation in halos formation 
seems neither to yield better constraints on the DM properties, nor to solve the slight tension 
between CMB and Gunn-Peterson data. We thus essentially confirm similar results obtained 
in the past, see e.g. [22], Planck data are not changing these conclusions in any significant 
way. On one hand, this reassures about the robustness of the reionization constraints to 
DM obtained by considering only the smooth contribution. On the other hand, this suggests 
that it is very hard to improve over them by including the relatively low-z contribution from 
halos. Barring very different particle physics or halo assembly histories (for some example see 
e.g. [49]), it appears that the only hope to revisit this conclusion in the future and to reveal 
some contribution of DM halos would be to measure x e or Tm as a function of z, especially at 
high redshift (z > 10), and at the same time, to improve our knowledge (both theoretically 
and observationally) on reionisation by the first stars. 

4 Impact of reionization histories on the CMB spectra 

In this section we go one step further in the discussion of CMB sensitivity to different reioniza¬ 
tion models, beyond the simple integral constraint on r re i 0 discussed in the previous section. 
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While unnecessary to settle the issue of current sensitivity to DM halo signals, this is useful 
to assess the capabilities to improve over current constraints with future CMB data. 

A first important point to make is that the CMB spectrum is in principle sensitive to the 
entire ionization history x e (z). The main effect of DM annihilation is to delay recombination, 
to increase the width of the last scattering surface, and to enhance the number of residual 
free electrons after decoupling. A larger ionization fraction x e results in more Thomson 
scattering of photons along the line-of-sight. As long as x e remains close to its asymptotic 
freeze-out value, rescattering impacts CMB observables at a very small level (although not 
totally negligible). When stars and/or DM halo formation start, x e increases with time, and 
the fraction of rescattered CMB photons becomes significant. Temperature and polarization 
anisotropies are damped on sub-Hubble scales, and regenerated on scales comparable to the 
Hubble radius. 

Typically, for usual models assuming single-step reionization in the range 6 < z < 12, 
high Z’s probe only the integrated parameter r re i 0 , since the temperature and polarization 
anisotropy spectra are mainly suppressed by a factor e^ 2Treio (this behavior can be spoiled to 
some extent by a more complicated reionization history, see appendix D for details). However, 
small Z’s are sensitive to the full reionization history - especially as far as polarization is 
concerned. Even in single-step models with a fixed r re i 0 , the shape of the CMB spectra in the 
range Z ~ 20 — 40 keeps an imprint of the details of the reionization history. In presence of 
DM annihilation in halos, x e tends to increase slowly at higher redshift, and a wider multipole 
range Z ~ 20 — 200 can in principle be impacted. 

In Fig. 6 we show the temperature (left) and EE polarization (right) multipole spectra 
computed for the two models of astrophysical reionization described in Sec. 2.1, with the 
bottom panels showing the relative difference between the two models, compared with the 
cosmic variance (shaded areas). They have been produced with CLASS version 2.4.3 for a 
few DM models with fixed 0 S = 1.04077 X 10“ 2 and r re i 0 = 0.079, in agreement with Planck 
measurement (TT,TE,EE+lowP [25]). The same 6 S is obtained by adjusting Hq (with fixed 
w&, Wcdm)- The fact that the variations shown in these plots are well within the limits of cosmic 



Figure 6. Upper pannels — Cj TT , C EE for the two models of astrophysical reionization de¬ 
scribed in Sec. 2.1. The shaded area represents cosmic variance. The baseline ACDM model is 
assumed, with T re j 0 = 0.079 and 9 S = 1.04077 x 10 -2 fixed in agreement with Planck measurement 
(TT,TE,EE+lowP [25]). Lower pannels — Relative difference between the two models. 
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variance exemplifies why, despite the fact that a sensitivity to x e (z) is present in principle, 
it is considered to be hopeless to infer information on different astrophysical reionization 
scenarios via CMB observations, and why to a large extent in this framework it is an excellent 
approximation to assume that CMB is only sensitive to r re i 0 . Nonetheless, note that a greater 
sensitivity to x e (z) of the EE polarization with respect to TT spectrum is still manifest. 




Figure 7. As in Fig. 6, but with the DM effects (both smooth and in halos) now added, for a single 
reference case of halo boost factor. 




Figure 8. Zoom on the residuals with respect to the baseline models reported in Fig. 6, with several 
assumptions for the halo contribution reported and for two choices of the astrophysical reionization. 

The situation is partially altered if DM annihilation in halos is added. In Fig. 7, we show 
plots similar to Fig. 6, with the DM effects (both smooth and in halos) now added; we fix the 
DM channel to the muon one, the DM mass to 1 GeV, and p ann = 2.3x 10~ 7 , saturating current 
bounds. While in Fig. 7 only a single, realistic value 10 6 for the boost factor is shown, in Fig. 8, 
we present the residuals with respect to the reference ACDM models, with growing effect of 
the halo term (besides the 10 6 one, also a factor 10 and 100 larger). It is clear that, at least for 
the phenomenological model for the astrophysical reionization, potentially detectable effects 
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emerge. Since we fixed 9 S and r re i 0 , we nearly eliminated any oscillatory patterns and step-like 
discrepancy between the set of curves describing the TT spectrum. This is not entirely true 
in models with a significant effect of annihilation in halos, for which a residual step-like effect 
can be clearly observed in Figure 8. This has to do with the ambiguity in defining r re i 0 in 
presence of partial reionisation at high-redshift, as explained in Appendix D. By just fixing 
Aeioj we do not eliminate completely the dependance of the high-/ temperature spectrum 
on the reionisation history. Apart from this effect, which could be compensated by a shift 
in the reionisation time or in the overall normalisation of the primordial spectrum, changes 
remain very modest and limited to small scales (where degeneracies with other parametrical 
extensions of ACDM are likely present). On the other hand, they can become sizable and 
characterized by a peculiar pattern at low-/ on the polarization power spectrum, coming from 
the enhanced rescattering of photons at intermediate redshifts z re i 0 < z < 600. While for 
realistic values of the halo parameters and DM mass it is still likely that this falls within the 
cosmic variance, it cannot be excluded that more extreme values of the parameters could be 
independently ruled out by CMB polarization observables, provided that forthcoming final 
release of Planck data manages to keep under control the systematics at low—/. 

5 Discussion and prospects 

The CMB temperature and polarization anisotropy pattern is sensitive to the energy injection 
by DM annihilation, essentially only via its constraints on the reionization optical depth r re i 0 
(we clarified in appendix D in which limit this is actually true). In turn, the optical depth r re i 0 
probed by the CMB depends on the reionization history, i.e. the function x e (z), which itself 
couples to the thermal history of the gas, controlled by the function T\/\{z) (their standard 
evolution equations have been reported in Appendix A for completeness). One major element 
of novelty of this work has been to study for the first time the dependence of the cosmological 
observables x e (z) and Tm(z) from the underlying astrophysical reionization/heating source 
model, both in presence and absence of DM sources. We modified the dynamical equations for 
x e and Tm (see Sec. 2.1) arguing that virtually all previous treatments have been incomplete 
and inconsistent in that respect. 

We have then revisited the problem of current CMB constraints on dark matter anni¬ 
hilation in halos, with a state-of-the-art treatment. We have followed the by now standard 
formalism of Ref. [13] and made use of the numerical tools provided in [19], to compute pre¬ 
cisely how the energy is deposited in the medium. We have clarified and corrected a few 
mistakes in previous work and improved over them. In appendix B, we provide a detailed 
comparison of the many formalisms used in the literature to study DM annihilation in halos 
and carefully explain where are the few mistakes or approximations over which we have im¬ 
proved. Appendix C describes in details the parametric model of DM halo formation adopted 
in our calculations, although a generalization to different parameterizations is straightforward 
and would not alter our conclusions. 

In agreement with most previous literature, we have confirmed that with conventional 
assumptions on DM models and halo formation, DM annihilation fails to play a dominant 
role in reionizing the universe, and can at most coexist with reionization from stars. Only 
very unrealistic halos could give a significant contribution to the reionization optical depth, 
and even then, such models are hardly compatible with astrophysical measurements of the 
ionization fraction at low redshift. No plausible DM model can produce a too high reionization 
optical depth, hence CMB measurements of r re i 0 do not provide additional constraints on the 
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annihilation cross section. At face value, this means that previous constrains derived assuming 
DM annihilation in the smooth background only are robust and independent of uncertainties 
on structure formation. Note that the recent update of the tools to compute the energy 
deposition in the medium provided by [24], following the results of Ref. [21] would not alter 
our conclusions: since the account for the new channel of energy loss through very low energy 
photons (with energy < 13.6 eV) produced during the development of the electromagnetic 
cascade, such a refinement would only make the impact of dark matter annihilations in halos 
even slightly weaker. We have checked this conclusion following the so-called "approximate" 
method described in Refs. [21, 24] which consists in withdrawing some power to the transfer 
functions of [19], using a specific new table supplied in Ref. [24], This method is, according to 
the authors, as precise as the new transfers functions and is better suited for our treatment. 

To assess if these rather pessimistic conclusions are linked to the intrinsic insensitivity 
of the CMB to these effects or merely to the current lack of precision, in section 4 we have 
computed the CMB angular power spectra, and compared the results with or without an¬ 
nihilation in halos. We have shown that, within standard assumptions for the astrophysical 
reionization and for plausible halo models, both the effects on the TT and EE multipole spec¬ 
tra are unobservable, falling below the level of cosmic variance. One would conclude that fits 
of CMB spectra in presence of DM annihilation in halos cannot provide (even in principle) 
better constraints on DM models. However, we have also shown that to some extent this 
conclusion can be altered if a different scenario for the astrophysical source of reionization is 
adopted. At present, it cannot be excluded that some extreme but viable DM halo parameter 
space might be eventually probed by CMB polarization data, provided errors can be kept at 
the level of the cosmic variance. To the best of our knowledge, it is the first time that this 
effect is highlighted. 



Figure 9. As in Fig. 4 (see also the right panel in Fig. 3) but for the IGM temperature. 
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Looking beyond the CMB probe, our results also leave the door open to some encouraging 
perspectives to further constrain DM annihilation in halos via other obserbables, more directly 
linked to x e (z) or Tm{z). To illustrate this point, one can look for instance at Fig. 9, where 
we report the evolution of Tm for several mixed reionization models, analogous to what 
shown in Fig 4 for the evolution of x e (z). In particular, green lines are benchmark evolutions 
for purely astrophysical scenarios (dashed: single step reionization/reheating; dot-dashed, 
phenomenological one; solid line: standard case where only the feedback from x e (z) is included 
in the dynamics of T m ). The blue and red versions of the corresponding lines show the case 
where both smooth DM injection and halo one have been added as well, with growing role 
of halos: DM annihilation in halos with z h = 20, B{z = 0) = 10 6 * are reported in blue, and 
Zh = 30, B(z = 0) = 10 8 in red, respectively. 

The dot-dashed red line shows that already with current rudimentary constraints, if one 
could trust the phenomenologically motivated model for the astrophysical heating, extreme 
halo parameters could be excluded. At the same time, as shown by the solid green curve, it is 
immediately obvious that current treatments of the Tm(z) evolution are strongly inadequate, 
and a significant effort should be put in achieving a realistic modeling of the sources of heating. 
The qualitatively most interesting effect of DM annihilation in halos is the possibility of a 
sign shift in the difference between T 7 and Tm at much higher redshift than expected for 
astrophysical models. This could be detectable for instance with 21cm surveys, since the signal 
would appear either in absorption or in emission (e.g. [50] and reference therein). Tomographic 
surveys of the cosmological 21cm observable, sensitive both to x e (z) and Tm(z), are certainly 
the most promising avenue to progress in the knowledge of these redshift epochs. Numerous 
experiments such as PAPER 64 6 , 21CMA ', MWA 8 , LOFAR 9 , HERA 10 or SKA 11 , are 
now (or will be) attempting at measuring the 21 cm signal. Hopefully, it could be possible to 
see the impact of halos through such a sign shift at relatively high redshift, when the stars are 
not yet “polluting” the signal. If not, a good description of the influence of stars on Tm and x e 
would be of first importance, as we have illustrated using semi-empirical models taken from 
[36, 38, 51, 52], Some authors are also trying to improve on the modelling of these effects 
using numerical simulation (such as 2 lcm FAST [53-55]). Achieving a sufficiently accurate 
treatment of the interplay of stellar and exotic sources is definitely a complicated task and a 
long term goal, but certainly it deserves further investigations. 

Acknowledgments 

We thank Marco Cirelli, Sergio Palomares Ruis and Aaron Vincent for very useful discussions. 
We also thank the anonymous referee for very constructive comments. Support by the Labex 
grant ENIGMASS is acknowledged. 

A Appendix A: Peebles’s recombination 

We want here to briefly introduce the standard tools useful for computing observables rel¬ 
evant to the recombination and reionization epoch. The fundamental equations have been 

6 http://eor.berkeley.edu 

' http: //21cma.bao.ac.cn 

8 http://www.mwatelescope.org 

9 http://www.lofar.org 

10 http://reionization.org 

“http://www.skatelescope.org 
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introduced by Peebles, quickly followed by Zel’dovich, Kurt and Sunyaev in the seminal pa¬ 
pers [56, 57]. Since then, there has been substantial improvements of the model, as detailed 
in [33, 34, 58, 59] and references therein. In the following, we will briefly sketch the basic 
ideas in order to keep our article self-contained. 

In the standard picture, neutral hydrogen is not directly formed by capturing an electron in 
the ground state (Is), because this process yields a photon of the exact energy necessary to 
reionize a neutral atom, and hence results in no net decrease of x e . For a similar reason the 2 p 
level gives no significant net capture, since the emitted photons will re-excite another atom, 
that can then easily escape if the atom interacts with a CMB photon. Recombination hence 
typically requires an electronic capture in the 2s excited state, with a de-excitation through 
a forbidden transition, in order for these photons not to re-ionize other atoms. The accurate 
calculation of this process requires to follow coupled kinetic equations for the free electrons 
and thermal photons (taking into account a certain number of hydrogen levels), a problem 
typically solved numerically. The most recent papers about recombination [60, 61] distin¬ 
guish between HII, Helll and Hell recombinations. They have highlighted very interesting 
effects on the CMB spectrum itself, through deviations from the black-body shape known as 
the “recombination lines”, which are in principle detectable with a dedicated experiment, but 
such refinements are not necessary for our purposes. 

Assuming that helium has entirely recombined, which is the case in the redshift range 
we are interested in, the evolution equation for the free electron fraction x e = n(HII)/n(H) = 
n(HeII)/n(He) in term of the redshift z writes: 

dx e (z) 1 /n/ , T , ,, . . 

~d^ " (1 - z)H(z) ^ 


with 


R{z) = C 


anxlnn , 


I{z) = C 


his g 

/ 3 f/(l — x e )e k b T M 
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These two terms are respectively the standard recombination and reionization rates. The 
first term encodes the probability that one free electron encounters an ionized hydrogen, is 
captured but not directly in the ground state and finally decays from the n = 2 state to 
the ground state before being ionized. The second term encodes the probability that a CMB 
photon redshifts at the Lyman-a transition frequency, and hits a neutral hydrogen in the 
n = 2 state. This standard scenario is known as Peebles "Case B" recombination. In this 
framework, the coefficient C is given by: 


1 + 7C g A g n g (l - x e ) 

1 + Kh{Rh + /3ff)n#( 1 - x e ) 


(A.3) 


where A// = 8.22458 s 1 is the decay rate of the metastable 2s level, n#(l — x e ) is the number 

A 3 

of H atoms in the ground state, and Kn = accounts for the cosmological redshifting 

of Lyman-a photons. Finally, an and fin are the effective recombination and photoionization 
rates for principal quantum numbers > 2 in Case B recombination (per atom in the 2s state), 
v a is the Lyman-a frequency and Tm is the temperature of the baryonic gas. 

The evolution of the matter temperature is instead ruled by the equation: 
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where the dimensionless parameter 7 , also called opacity of the gas, is defined as : 

_ 8(TTarTQ MB x e 

3 Hm e c 1 + f He + x e 

with <tt the Thomson cross-section, a r the radiation constant, m e the electron mass, c the 
speed of light, and /# e (= 1^,/[4(1 — Y/)]) the fraction of helium by number of nuclei. The 
first term is an adiabatic cooling due to the Universe expansion (for a massive non-relativistic 
particle the kinetic energy E\. = p 2 /2m oc a -2 and can be related to the averaged energy 
of indivual particles in a monoatomic gas Ek = 3/2& 7 T), whereas the second one accounts 
for interactions between CMB photons and matter, mainly through Compton scattering with 
free electrons. The evolution of Tm therefore mainly depends on the second term, with the 
sign of the difference Tcmb — Tvi determining if collisions heat or cool down the gas, and the 
value of 7 responsible for the coupling strength: If 7 S> 1, the second term dominates and 
Compton scatterings thermally couple CMB photons and baryons (Tm oc Tcmb oc (1 + z)), 
whereas Tm adiabatically decays as Tm oc (1 + z) 2 when 7 <C 1 . 


B Appendix B: Remarks on the efficiency function /(z) 

The computation of the /(z) for each DM model (mass and final state) is the key ingredient 
needed to correctly describe impact of DM annihilation on the cosmic gas. For typical WIMP 
models, this computation was made in [13, 19, 21], and here we will only sketch the procedure. 
We refer reader to previous references for detailed explanations and estimates of the systematic 
uncertainties relative to the treatment. 

In order to compute these functions, one needs to follow the evolution of an electromagnetic 
cascade starting at the DM mass (typically O(GeV-TeV)) until the energy at which atomic 
processes dominate (sub-keV) and over the entire redshift range that we are interested in, 
typically from z ~ few thousands until today. This makes > 10 orders of magnitudes in 
energy and 3 orders of magnitude in redshift that the code should be able to span. 

A significant simplification can be achieved by splitting the cascade in two regimes: at 
high energy, interactions with CMB photons dominate and are responsible for particle mul¬ 
tiplication and energy degradation; at low energies, electromagnetic particles interact with 
the gas and lead to heating, excitation and ionization of the medium. In practice, at each 
redshift, the first part of the computation provides the fraction of energy available to interact 
with the matter (which, going beyond the on-the-spot approximation, can be bigger than 
1 ), namely the efficiency factor /(z), whereas the second part of the computation yields the 
energy repartition fractions %j- 


The last remark we want to make consists in a comparison of the formalism used by 
authors to implement the energy deposition from DM annihilations. Instead of parametrizing 
the energy deposition history through the /(z) functions, the authors of [10, 11, 14, 20] have 
used a more explicit method, which has however some caveats. Following notations of [11], 
equivalent to other articles, one can write after a few manipulations the energy deposition as: 
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where (t(E^) is obtained by summing cr lA _^ e - A+ , cr^ e -^ e - and a lA _^ e ± A (the index A refers 
to H and He atoms), n A (z) scales as (1 + z) 3 , (^7( El { ) is the spectrum of photons produced by 
the DM annihilations (both prompt and from inverse Compton scattering), El. is the energy 
at the time of injection and E 7 is the energy at the time of deposition, and k(z' , z; El f ) plays 
the role of an absorption coefficient defined as: 

C z rlf, 

k(z ,*',£ 7 )~jdz"—n A {z")a{E”) . (B.2) 

Three comments can be made about this formula. 

• First of all, a mistake is made in eq. (4) of Ref. [14] (the equivalent of our eq. (B.l)) and 
that mistake propagated to Ref. [20]. Indeed, in eq. (4) of Ref. [14], the last integral 
is performed over E 1 (the energy at the time of injection) instead of E 7 (the energy at 
the time of deposition). Since the two are related by E' = Ey( 1 + z')/( 1 + z), the final 
result differs by a factor (1 + z) 2 / (1 + z') 2 in the last integral over dz'. This can be seen 
explicitely in eq. (4.21) of Ref. [20], where the exponent 6 should in fact be 8 . 

• Secondly, in order to perform the integral analytically, Refs. [14, 20] assume that it is 
possible to consider that all interactions between injected photons and matter are due 
to Compton scattering in the Thomson limit. This approximation works well for the 
smooth background, but it is inaccurate compared to the approach of Slatyer et al. [13] 
in the context of annihilation in halos. 


• To further compare the approach of Refs. [10, 11, 14, 20] with that of Ref. [13], it is 
possible to manipulate eq. (B.l) in such way to absorb all the energy and redshift de- 
pendance. Then, the equivalent of the energy deposition function f(z) appears, defined 
as: 


dE 

dVdt 


/ \ 7 / \ 2/ \ wann'W, 

(z) = f{z)p (z)- 

dep m DM 
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where 


/(*) = 


s: 


dz'(l + z) 5 n A (z)p 2 (z') 
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Even if this formalism goes one step beyond the “on-the-spot” approximation thanks to 
the integral over z\ the complexity of the cascade evolution is not taken into account. 
The two approaches do not give significantly different results when considering only the 
smooth background, since the CMB experiments are not sensitive to the shape of f(z) 
at low z, but using eq. (B.4) can lead to wrong results when halos are turned on. 


C Appendix C: The boost function B(z) 


Generally, it is possible to re-write the energy deposition rate in halos by introducing an 
enhanced dark matter density p 2 (z) = (1 + I3(z))p 2 (z). One now has: 


dE 


2 2 
= Pc c 


^dm (1 + z ) 
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(1 + B(z)) 


dVdt 
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To compute this "boost" of the density, several ways have been introduced. In the so-called 
"Halo model" (HM) framework, it is assumed that all the mass in the Universe is contained 
in virialized objects (the halos), fully characterized by their mass. The key ingredients are 
the spatial distribution of matter inside a halo (usually called the density profile ) and the 
number of halos per unit mass (namely, the mass function, that can evolve with the redshift). 
This allows one to express the boost B(z) as : 


m= i¥hu ii+z) (C2) 

where M m ; n is the minimal mass of DM halos and r2oo the radius of a sphere enclosing 
a mean density equal to 200 times the background density (some authors prefer to use the 
virial radius r v , but they are strictly equivalent). The last integral is usually recast in terms of 
the concentration function F(ch), which depends on the concentration parameter Ch = V2oo/ r s 
with r s the scale radius of the given profile: 



drAirr 2 p\(r) 


Mp(zF ) f 

3 \ n M 


2 

F(c h) 


(C.3) 


where zp is the redshift of halo formation (not yet well known) and p(zp ) = 200 p c Hm(1 + ~f) 3 
the average matter density within a sphere of radius r 2 oo- 

One now needs to relate the concentration function to the halo profile, and several types of 
profile have been proposed in the litterature. For instance, the most commonly used are the 
Navarro-Frenk-White (NFW) profile [62], Einasto profile [63, 64] and isothermal-like Burkert 
profile [65], the former two showing a more peaked distribution than the latter. The last key 
quantity - the mass function - which in this framework also encodes how the density evolves 
with the redshift, can be calculated analytically in the Press-Schechter formalism [66], leading 
to : 


dn 

dM 


p c Q m d In a 1 
dM 


fi?) 


(C.4) 


where the variance of the density held a and the differential mass function /(<?) have been 
introduced. We shall not develop further the HM and the Press-Schechter formalism and 
refer to [20] for more details. Following their development, one can show that, in the frame 
of the HM with the Press-Schechter prescription for the mass function, the boost from halos 
is given by: 


B(z) 


fh 

( 1 + z ) 3 


erfc 


1 + z\ 

1 + ZhJ 


(C.5) 


where fh = ^(1 + ZF) 3 F(ch) and Zh is the characteristic redshift at which halos start 
to contribute (we typically have Zf — 2%). It could be possible to use more advanced 
prescription for the halo mass function, such as the Sheth-Tormen formula [67] or the one 
from Tinker et al. [68], but this would not lead to major differences for the problem we are 
dealing with. Indeed, as it appears in this framework, the mass function mainly dictates 
the evolution of the distribution with z and is only known for very small redshifts compared 
to the range we want to span. The error function is merely a prescription inspired by a 
simple model to describe the transition from a smooth distribution of matter to structures 
of virialized objects, whose normalization at low redshift is encoded in the function ff. To 
learn how structure formation precisely happens, especially at high z, is one ingredient that 
would allow one to improve over the present study. Here, we will limit ourselves to vary 
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z\y in the interval [20,30] to ease the comparison with existing literature. Concerning the 
additional parameter /(,, we deduce a reasonable range for it from the range of B(z = 0) 
computed in [47, 48]. In those references, an alternative method is used to compute the boost 
factor, by direct integration of the power-spectrum down to very small scales (where a cutoff 
is induced by WIMP dark matter free-streaming or kinetic decoupling) calibrated directly 
on simulations. The large uncertainties in dealing with the power-spectrum in the deeply 
non-linear range explain the broad interval B(z = 0) 6 [10 4 ,10 8 ] inferred for typical WIMP 
candidates, although more recent simulations tend to narrow it down to [5 x 10 4 ,10 6 ]. 


D Appendix D: Discussion on r reio as it is measured by Planck 


The Planck collaboration explains in sec 3.4 of Ref. [25] how they proceed to obtain a mea¬ 
surement of the optical depth to reionization r re i 0 - Let us recall briefly how this is done, and 
why this might be problematic when one modifies the reionization history in a non-trivial 
way and wants to compare the new r re i 0 to the Planck results. 

In general, the optical depth at redhift z is defined as 


r(z) 


f z rH- 

I n H {z)x e (z)a T —dz' . 


(D.l) 


Planck obtains an estimate of the reionization optical depth by: (i) assuming single- 
step reionization where x e (z ) is described by a postulated function centered on an adjustable 
redshift z re i 0 ; (h) computing the corresponding CMB temperature and polarization power 
spectra using the Bolztmann codes CAMB [43] (or CLASS [32, 69] for cross-checks); (iii) 
comparing these spectra with those extracted from observed maps. The fit directly gives 
some bounds on the reionization redshift z re i 0 - At the same time, for each model, the codes 
also compute the integrated parameter T^™ b = r(z c ), where the cut-off redshift z c takes an 
arbitrary default value z c = 40 in both CAMB and CLASS, and reports bounds on this T^[™ b . 
Actually, since there is a one-to-one correspondance between T^[™ b and z re io (valid for the 
postulated category of reionization models), the Bayesian parameter extraction can even be 
done with a flat prior on rather than z re io- 

In such models, reionization affects: 


• the low-/ part of the polarization spectra, creating the so-call reionization bump, 

• the high-/ part of the temperature and polarization spectra, which are step-like sup¬ 
pressed, with a suppression factor saturating above / ~ 70 at a value very close to 

_9_camb 

6 Z reio 


The first impact is very challenging to detect because of the smallness of the signal, which 
requires a very good control of instrumental systematics and polarized foreground emissions. 
Even if theoretically, the signature of reionization on the CMB polarization spectra is very 
clear, the actual measurement of the El-mode low-/ spectrum still has quite large uncertainties, 
and that of the R-mode spectrum is not sensitive enough to probe any reionization bump. 
Hence, current sensitivity to reionization comes mainly from the second of the two effects. 
The data probe mainly the integrated parameter T r c J]™ b , and the bounds reported on this 
parameter are often thought to be model-indepent (i.e., valid for any reionization history). 

However, some of the above statements are only true for models such that the ion¬ 
ization fraction x e (z ) only starts raising at low redshift. For more complicated models like 
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those involving DM annihilation, the fact that there is a step-like suppression of the C/’s, 
and that the suppression factor asymptotes to e _2Trei ° , are not necessarily accurate. For 
instance, Figure 10 shows that in presence of DM annihilation in halos, the reionization effect 
is not exactly step-like. It has superimposed oscillations, and the suppression factor drifts 
significantly between Z ~ 100 and Z ~ 2500. 



multipole l z 


Figure 10. Left pannel - Comparison of the CMB temperature power spectrum for several mixed 
star-DM reionization models, with z re ; 0 = 6.5 and a DM particle of 1 GeV annihilating into electrons, 
normalized to the CMB power spectrum in a universe without reionization. Right pannel - The 
corresponding ionization histories. 

This can easily be understood using the line-of-sight integral formalism [70], which shows 
that the Sachs-Wolfe and Dopper terms in the CMB transfer functions are given by an in¬ 
tegral over time of a product of cosmological perturbations, Bessel functions, and finally the 
visibility function g = t' e~ T (as well as its time derivative in the Doppler term). The visibility 
function g(z) gives the probability of last scattering, and is normalised to one by construc¬ 
tion: f g(z) dz = 1. In absence of reionization, g{z) has a single peak at recombination, and 
it could be approximated by zero for £ <C Zdec without changing the result (physically, this 
means that all observed photons last scattered near z = z^ec)- With conventional single-step 
reionization model (black line in fig. (11)), a second peak centered near z = z re i 0 is visible. 
For high Z’s, the integral over time still only picks up contributions near z = Zd ec - However, 
the normalization of g(z) implies that the amplitude of the recombination peak in g(z) is 
reduced by exactly e~ T ^ z \ where z can be choosen anywhere in the range z re i 0 <C z <C z ^ ec 
(e.g. z = 40, as taken by default in CAMB and CLASS). Hence the high-Z temperature and 
polarization power spectra are supressed by exactly e _2T reio . 

This conclusion is almost unspoiled by the introduction of DM annihilation in the smooth 
background, or for relatively small halo contribution (blue and green curves in fig. (11)). 
However, the non-standard shape of g in the case where reionization in DM halos is very 
strong and starts early (cyan and red curves in fig. (11)) produces non-trivial effects, because 
the line-of-sight integral starts to pick up sizeable contributions coming from the second 
reionization bump. Physically, this comes from extra small-scale anisotropies generated by 
CMB photons rescattering at medium redshift, for instance between z ~ 40 and z ~ 100. 

Since Planck measures the optical depth mainly throught its suppression effects on the 
C{ s at high Z, it is possible to say that what Planck really measures is an effective parameter 
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Figure 11. Visibility function for the reionization history from fig. 10, compared to the visibility 
function in the absence of reionization. 


r e fj = — ^ log {Cf T /Cj T ' no rci °), where l is some effective multipole value in the region where 
the data has maximum sensitivity, and we are assuming that most of the sensitivity comes 
from temperature. For usual single-step reionization models, r e fj coincides both with T^™ b , 
computed up to z c = 40, and with a more sensible and robust definition that one could give 
of the reionization optical depth, where the integral would run up to the redshift at which 
the ionization fraction is minimal: 


T reio =t(z*) , *'(*•) = 0 . (D.2) 

For models with a raise in x e (z) at high redshift, it is obvious that r e ff and r r C g™ b can be very 
different, because integrating up to z c = 40 is not sufficient. We also checked this numerically 
and found a significant mismatch. A more interesting question is whether r e ff and r™" are 
close to each other. In other words, can we still say that the total effect on the CMB spectra 
at high l is related to the optical depth computed up to the very beginning of reionization? 
We tested this Ansatz explicitly, and found that the difference between these two quantities 
is small, being of the order of 10% in the worst cases, i.e., still within the error bars of Planck. 
In summary: 


• for usual single-step reionization: 


camb 
' reio 


7"eff 7"n 


• in more general cases with an early enhancement of the ionization fraction: 


camb 

'reio 




_ _ . ^min 

Teff — Aeio ■ 


Therefore, as a first approximation, it is safe to define everywhere in our analysis r re i 0 as r™ ™. 
and to compare it with the value of r re i 0 reported by Planck. This conclusion legitimates the 
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analysis performed in section 3. For even better accuracy, the next step would be to perform a 
full parameter extraction in which exact power spectra are fitted to the data, and to derive new 
bounds on r re i 0 valid in all cases, but this is not necessary for reaching the main conclusions 
of this paper. 
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